Rectification of Vortex Motion in a Circular Ratchet Channel 
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We study the dynamics of vortices in an asymmetric (i.e., consisting of triangular cells) ring 
channel driven by an external current / in a Corbino setup. The asymmetric potential can rectify 
the motion of vortices and cause a net flow without any unbiased external drive, which is called 
ratchet effect. With an applied ac current, the potential can rectify the motion of vortices in 
the channel and induce a dc net flow. We show that the net flow of vortices in such a system 
strongly depends on vortex density and frequency of the driving force/current. Depending on the 
density, we distinguish a "single-vortex" rectification regime (for low density, when each vortex is 
rectified individually) determined by the potential-energy landscape inside each cell of the channel 
(i.e., "hard" and "easy" directions of motion) and "multi- vortex" , or "collective", rectification (high 
density case) when the interaction between vortices becomes important. The frequency of the driving 
ac current determines a possible distance that a vortex could move during one period. For high 
frequency current, vortices only oscillate in the triangular cell. For low frequency, the vortex angular 
velocity uj increases nearly linearly until the driving force reaches the maximum friction force in the 
hard direction. Furthermore, the commensurability between the number of vortices and the number 
of cells results in a stepwise u — I curve, which means that the average angular velocity cj of the 
vortices is a discontinuous function of the applied current. Besides the "integer" steps, i.e., the large 
steps found in the single vortex case, we also found fractional steps corresponding to fractional ratio 
between the numbers of vortices and triangular cells. The principal and fractional frequencies for 
different currents are found, when the net flow of vortices reaches the maximum that is proportional 
to the frequency when the density of vortices is low. We have performed preliminary measurements 
on a device containing a single weak-pinning circular ratchet channel in a Corbino geometry and 
observed a substantial asymmetric vortex response. 

PACS numbers: 74.25. Uv, 74.78.Na, 05.70.Ln 



I. INTRODUCTION 



A net flow of particles under unbiased external fluc- 
tuations/drive due to an asymmetric potential, which is 
called ratchet effect, has received much attention during 
the last decades. The transport and dynamical prop- 
erties of particles on asymmetric potential have been 
widely studied, e.g., in physics and biology [iHgl- Ran- 
dom motion of particles can be rectified in such an asym- 
metric system, which can be used for, e.g., controlling 
particles motion, separating different types of particles 
(i.e., molecular sieves), for both underdamped and over- 
damped particles Q and for molecular motors [5|. Vor- 
tices in a type II superconductor often (e.g., for mag- 
netic field close to Hd) can be treated as classical over- 
damped "particles" . Most of the experiments on vortex 
motion rectification used arrays of asymmetric pinning 
sites (e.g., nanoengineercd antidots or triangular mag- 
netic dots/inpurtics) to create an asymmetric potential, 
which rectifies the motion of vortices [8l-[l4l. |27| . The rec- 
tified vortex motion was directly observed in experiments 
by imaging vortices via Lorentz microscopy [l5|. Peri- 
odic arrangement of point defects of a gradual density or 
periodic square array of ferromagnetic dots of decreasing 
size, i.e., varying the density of pinning sites or the size of 
pin ning sites, were shown to result in a ratchet potential 
[la, ll3- When vortices are trapped by pinning sites, the 



repulsive vortex-vortex interaction creates a higher en- 
ergy barrier near the area with higher density of pinning 
sites. Therefore, an asymmetric potential can be created 
by the gradient of tlie density of pinning sites |16| . Even 
without spatial asymmetry (i.e., without any asymmetric 
walls/boundaries or asymmetric pinning sites), the mo- 
tion of vortices still can be controlled by time-asymmetric 
driving force p8l - l20l |. Due to the possibility of control- 
ling their motion, the dynamical behavior of vortices in 
such systems has attracted considerable interest. A series 
of elastic and plastic vortex flow phases were found |2l| - 
[23j . Besides the liquid- like and solid lattice phase, vortex 
motion also revealed a jamming behavior [2j, [2J]. When 
the density of vortices is changed, the vortex flowing di- 
rection can change to the opposite 0, [lO, ll3: l25| , which 
means vortices can drift in either the "hard" direction 
or the "easy" direction of the ratchet, depending on the 
vortex density. By controlling the motion of vortices, it is 
possible to remove vortices or reduce the vortex density 
by using a combination of two opposite oriented ratchet 
arrays [26[ . The order of vortices and commensurability 
between vortices and cells also play an important role in 
vortex dynamics [13, llj, l25|, [23 . In two-dimensional (2D) 
ratchets, the dynamics in the transverse ratchet was also 
studied in theory [4| and in experiments [23, [231 ■ 

In the present paper, we study the dynamics of vor- 
tices in a circular channel formed by asymmetric trian- 



gular (funnel) cells (TCs) (see Fig. [T]) [note that earlier 
this approach to form asymmetric channels in experi- 
ment (i.e., using weak-pinning channels) was employed 
in a stripe geometry [6|]. Due to the radially flowing 
current in a Corbino setup, the driving force is not uni- 
form inside the cell which is different from linear ratchet 
channels. Such a geometry (i.e., asymmetric in the az- 
imuthal direction and in the radial direction), as will be 
shown, leads to a specific dynamical behavior, for ex- 
ample, a vortex located near the inner corner of a TC 
(which is closer to the center), experiences a stronger 
driving force and moves to the next TC while a vortex 
located near the outer corner of TC does not move. The 
circular geometry of the ratchet channel is convenient for 
studying commensurability and step-motor (phase lock- 
ing) behavior [IJ, [2J]. We analyze in detail low and 
high density regimes of rectification, i.e., "single- vortex" 
and "multi-vortex" regimes. We demonstrate that the 
mechanism of rectification are qualitatively different for 
these two cases. In addition, we have performed prelim- 
inary measurements using a single nanofabricated weak- 
pinning ratchet channel of a-NbCc with strong-pinning 
NbN channel edges in a Corbino set-up, and we observed 
a substantial asymmetric vortex response. 

This paper is organized as follows. First in Sec. II, 
the model of our systems will be presented. Then in 
Sec. Ill, we will show how the density of vortices and 
the frequency of driving current influences the dynam- 
ical behavior of vortices. Commensurability effects of 
vortex density and of the frequency of current will be 
discussed in Sec. IV. We present the measurements of a 
vortex ratchet in a Corbino geometry in Sec. V. Finally, 
the conclusions will be given in Sec. VI. 



II. MODEL AND SIMULATION 

We consider a ring-like weak-pinning channel con- 
structed of N partially overlapping TCs as shown in Fig. 
[TJ In our simulations, the radius of the ring R is typi- 
cally set as 6A, where A is the magnetic field penetration 
depth, the wider part of the channel (i.e., the base of 
TCs) w = 0.75A, and the width of the narrow part A 
(the neck) is typically 0.1 of the wider part (i.e., the ra- 
tio g = A/w = 0.1). We also performed simulations in a 
channel with a wider neck part, e.g., g = 0.15, for com- 
parison. The weak-pinning channel (where vortices can 
move freely) is surrounded by a strong-pinning super- 
conducting material [30| which is modeled by a medium 
where vortices cannot move. An external current / radi- 
ally flows from the center of the disk to the edge, result- 
ing in the density of current J{p) ^ I / p. Therefore, the 
closer a vortex to the center of the disk, the stronger the 
Lorentz force that acts on it. The driving force due to 
the radial current is (Sll-ISal 



where is the unit vector in the azimuthal direction in 
the disk plane, /o = $§/(27r/ioA'^) is the unit of force, 
-^0 = fJ-o^^I/i^od) is the dimensionless driving current, 
and $0 is the flux quantum. We perform Langevin-type 
molecular dynamics (MD) simulations and numerically 
integrate the overdamped equations of motion |35l438l | : 



di'i 



with 



f^-J2f^I + f^ + f' + f' 



(2) 



(3) 



where 77 is the dimensionless viscosity coefficient which is 
set here to unity. We note that in practice the viscosity 
coefficient [39| varies between different superconductors. 
As an example, parameters for a-NbGe, a typical weak- 
pinning material, correspond to 77 « 10~^ Ns/m^ [4G |. 
Using this value of 77 in our calculations results in typical 
maximum values of vortex linear velocity tj w 10^ m/s 
(for a l/im-thick film), which approaches, but does not 
exceed, the Larkin-Ovchinnikov critical velocity for this 
material H^-IH. 

The vortex interaction /^" is described by a first or- 
der of the modified Bessel function of the second kind 
Ki{rij/X), and the thermal force obeys 



and 



(flit)) = 0, (4) 



{fnt)fJit'))=2r^kBT6,,6it~t'). (5) 
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To model the vortex-boundary interaction, we assume 
an infinite potential wall at the boundary (i.e., vortices 
cannot leave the channel) which decays inside the trian- 
gular cell with the same dependence on position as the 
vortex- vortex interaction potential. The total interaction 
of a vortex with the channel boundaries is calculated by 
integrating the vortex-wall interaction potential over the 
geometrical boundary of the channel. The resulting po- 
tential due to the boundary is shown in Fig. [51 and the 
vortex-boundary force /j is directly calculated from the 
potential. 

In our simulations, we first set T > 0, when no current 
is applied, and then gradually decrease temperature to 
let the system to relax to the ground state. Then we set 
T = and apply an external driving, i.e., an ac current 
resulting in oscillating Lorentz force with frequency v and 
amplitude Jq that acts on vortices, to study the dynamics 
of the system. 



III. RECTIFICATION OF VORTEX MOTION 

A. Density of vortices 

The vortex density, i.e., the number of vortices per TC, 
plays an important role in the dynamics of the system. 




FIG. 1: (Color online) The geometry of the system. The 
widths of the wide part oj and narrow part A are shown in 
the figure, where g — A/w is the ratio between the two parts. 
Here g — O.f, the radius of the channel R = 6X and the 
channel is constructed by 36 TCs, i.e., A'^ — 36. 
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FIG. 2: (Color online) The contour map of the modeled po- 
tential. The lowest potential in a TC is close to the geometry 
center of the triangle and the potential near the boundary is 
high enough to prevent vortices from escaping the TC. 



For very low vortex density, the vortex- vortex interaction 
is small as compared to the interaction with the bound- 
ary, and therefore it can be neglected. With increasing 
the vortex density, the interaction between vortices be- 
comes important. Therefore, we roughly distinguish two 
regimes, i.e., a single- vortex regime and a multi- vortex 
regime, respectively, corresponding to low and high vor- 
tex density. 

In order to characterize the dynamical behavior of the 
system of vortices, we calculate the net angular veloc- 
ity of each vortex and take average over all the vortices 
for, e.g., 100 ac periods. The resulting average angular 
velocity (called further "angular velocity" ) is denoted as 
uj which is analyzed for different parameters of the sys- 
tem and drivings. We note that angular velocity can be 
related in a straightforward way to the flux-flow voltage 
[so, |31[ that one would typically measure in an experi- 
ment. 



One system contains L vortices and N triangles in a 
circular chain, e.g., in our simulations we take N = 36. 
If there is less than one vortex per cell, i.e., i < 36 in our 
system (see Fig. [T|), then the density of vortices is low 
enough and one can neglect the interaction between vor- 
tices. One can imagine that all the vortices are far away 
from each other and weakly interact with other vortices 
but they are strongly influenced by the ratchet poten- 
tial induced by the boundary. Therefore, the dynami- 
cal behavior in the low density case is similar to that of 
the system with just one vortex [shown in Fig. [3] (a)], 
which is considered in the single-vortex regime. Thus, 
for L — 1 — 36, the maximum of the angular velocity 
is the same value [see Fig. [3] (a)]. When the density 
increases and therefore the interaction between the vor- 
tices becomes appreciable, the maximum of the angular 
velocity utm first decreases to 0.05 for L ~ 55 [see Fig. [3] 
(b)] and then starts to increase [see Fig. [3] (c)] when the 
system turns to the regime of high vortex density. 

To understand this non-monotonic behavior of the 
maximum of the angular velocity a;,„ [shown in Fig. 
ini[d)], we studied the trajectories of moving vortices. 
When the system is in the single-vortex regime, where 
i < 36, all the vortices move along a nearly circular tra- 
jectory no matter whether the applied current is small 
or large, which can be considered as a one-dimensional 
(ID) motion (shown in Fig. \^. For low driving currents, 
each vortex oscillates near its initial position inside a TC 
[e.g., see Fig. 2] (a)], and when the driving force reaches 
some critical value, all the vortices move with a net an- 
gular velocity w in the easy direction [e.g., see Fig. |3] 
(b)]. However, in the multi-vortex regime, i.e., the high 
density case (e.g., L = 80), the motion of vortices is not 
ID any more. The ac current drives vortices to pass 
through the narrow part (i.e., neck) of the channel in the 
easy direction and when current is alternated, some of 
vortices are forced to move into the corner by the oth- 
ers and "freeze" thus blocking the motion in the hard 
direction [shown in Fig. [5]. Due to the asymmetry in 
the radial direction and radially decreasing current den- 
sity in the Corbino setup, a vortex near the inner corner, 
which is closer to the center of the disk, moves faster 
than the one near the outer corner. Therefore, at a spe- 
cific value of the current [e.g., at Iq = 9, as shown in 
Fig. [5ja)], the vortex near the inner corner has a larger 
angular velocity. When the current increases, more and 
more vortices move along the circular trajectory and the 
distribution of vortices becomes strongly inhomogeneous 
[shown in Fig. EJb)]. With further increasing current, 
this 2D trajectory becomes narrow and all the vortices 
move along the circle, i.e., the 2D motion becomes ID 
[shown in Fig. [5fc)] . With increasing density, more and 
more vortices remain in the area near the corners of TCs 
and block other vortices moving in the hard direction. 
Therefore, the maximum of vortex angular velocity in- 
creases. Even in the descending-velocity region of the 
single- vortex regime [see e.g., when Iq > 21 in Fig. [3Jc)], 
the angular velocity is still large. 
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FIG. 3: (Color online) The cj — Jo curves for different density 
of vortices: (a) the density of vortices is low, L = 1 — 36 (in 
the single- vortex regime), the angular velocity reaches the 
same maximum; (b) the density increases, L = 40 — 55, the 
maximum of the angular velocity decreases; (c) for further 
increasing density, i.e., high density of vortices, L = 55 — 90, 
the maximum starts to increase and the critical value of the 
driving current, when the angular velocity starts decreasing, 
becomes larger; (d) the maximum of the angular velocity ijjm. 
first remains the same until the number of vortices L > 36. 
Then the maximum ojm first decreases, but for L > 55 it 
starts to increase. 




FIG. 4: (Color online) Trajectories of vortex motion, in the 
case of the low density of vortices (L = 10) (a), (b), and in 
the intermediate case (L = 40) (c), (d), for different values of 
driving current: Jo = 6.5 (a), 13 (b), 9 (c), and 18 (d). The 
initial positions of the vortices are marked by gray circles and 
the final positions after one period of ac current are marked by 
black squares. The trajectories of vortices, which are plotted 
by solid lines, are nearly circular no matter that the current 
is large or small. 



The maximum of the angular velocity uj^ for different 
vorticities L is shown in Fig. |3l^d). In the single- vortex 
regime, where there is less than one vortex per TC, Wm re- 
mains almost the same value because of the weak interac- 
tion between vortices. When L becomes larger, there are 
more than one vortex per TC and the repulsive vortex- 



vortex interaction in the same TC makes the motion in 
either easy or hard directions much easier. Therefore, 
iOra first decreases when L > 36. However, for L > 55, 
the maximum tUm starts to increase. The reason for this 
behavior is the transition from ID motion to 2D motion. 
When vortices move along 2D trajectories (see Fig. [3]), 



the repulsive interaction between vortices still "helps" a 
vortex to move to the next TC in the easy direction but 
blocks the vortex motion in the hard direction until ap- 
plying a large enough driving current. 

The radial asymmetry facilitates the motion of vortices 
from inner/outer corners (driven by stronger/weaker 
force) one by one. This is different from the case of a 
linear ratchet channel where two vortices in the corners 
are driven by the same force and, as a result, arrive simul- 
taneously at the neck region leading to jamming. In an 
asymmetric channel, the symmetry is broken, and jam- 
ming occurs only for rather large driving force. 



that in the low/intermediate frequency case. Each vortex 
is localized in a specific TC, i.e., the vortex is only influ- 
enced by a single potential well but not by the ratchet 
potential. 

For comparison, we also calculate cj — /q curves for the 
case when a TC has a wider connection part {g = 0.15) 
[shown in Fig. IH^b), (d), (f)]. The function a;(/o) in 
general shows a similar behavior as for g — 0.1. However, 
for the intermediate frequency {v = 1) of the applied 
current, we obtain more steps in w — /q curve in the case 
of L = 1 [shown in Fig. [IJd)] than that for g = 0.1 [shown 
in Fig. m^c)]. This relates to the commensurability effect 
that will also be discussed in Sec. IV. 



B. Frequency dependence 

The dynamical behavior of vortices in the considered 
ratchet system is also strongly influenced by frequency 
of the ac current. In this subsection we study frequency 
dependence of the rectified vortex motion. For different 
frequencies i^, the w — /q curves are shown in Fig. [5] 

When the frequency is low, the dynamical behavior is 
similar for different vortex densities. As the drive ampli- 
tude /q is increased, the angular velocity first increases, 
reaches a maximum, and then decreases to zero. The 
curves, especially for L = 1, in Fig. |6l^a) are similar to 
the analytic result in a ratchet potential (e.g., see Fig. 
2 in Ref. [2g|)- I^i the case of L = 1, the first critical 
value of current Id [shown in Fig. |6l(a)] corresponds to 
the maximum friction in the easy direction /+ and the 
second one /c2 corresponds to the case when the driving 
force reaches the maximum friction in the hard direction 
/~. When L = 40, the first critical value of current de- 
creases to I'd, which means vortices are easier to move 
in the easy direction. If the density increases further, I'^^ 
decreases and 7^2 increases, for L = 80 [shown in Fig. 
[6lja)]. Therefore, the interaction between vortices allows 
vortices to move even easier in the easy direction and 
harder in the hard direction. 

For an intermediate frequency [e.g., i^ = 1 as shown in 
Fig. ini^c)], when the distance a vortex moves during one 
period is comparable to the size of a TC, the dynamical 
behavior of vortices will depend not only on the driving 
force but also on the vortex density. This is explained by 
the fact that the confinement force due to the boundary 
and the vortex- vortex interaction are also comparable. In 
this case we obtain different dynamical behavior in single- 
vortex and multi- vortex regimes which were discussed in 
the previous subsection. Due to the commensurability 
between the numbers of TCs and vortices, a step struc- 
ture of the u) — If) curve is revealed in the low density 
case (e.g., L = \ and L — 10), which will be discussed in 
more detail in Sec. IV. 

As one can predict, the angular velocity becomes zero 
[see Fig. |6l[e)] when the frequency is high enough be- 
cause a vortex oscillates near its initial position inside a 
cell and thus its motion cannot be rectified. Therefore, 
the variation of the angular velocity is much smaller than 



IV. COMMENSURABILITY EFFECT 

A. Commensurability of vortex density 

As we defined above, the system contains L vortices 
and N triangles in a circular chain. If there is a common 
integer (except one) in terms of which two numbers L and 
N can both be measured, then they are commensurate. 
Otherwise, they are incommensurate. Fig. [7|) shows the 
average angular velocity w as a function of /q for different 
commensurate ratios. The ascending part of the w(/o) 
curves (i.e., where w versus /q increases) is stepwise. Be- 
sides the large steps of angular velocity in the co — Iq curve 
for L = 1, which wc refer to as "integer steps", we also 
found smaller steps for some specific vorticities L (shown 
in Fig. [7]). If L/N — k/m and fc 7^ 1, where k and m are 
incommensurate integers, the small steps can be found in 
the uj — Iq curves [in Fig. [Il[b)-(d)]. When k = 1, we ob- 
serve only integer steps of the angular velocity a;(/o) [in 
Fig. EJa)]. The difference in the angular velocity between 
two adjacent steps is always ojq. If /c 7^ 1, then we can 
find a fractional step whose magnitude is cuo/k. For ex- 
ample, for fc = 2 [shown in Fig. [T^b)], the smallest step is 
ujo/2 in the system with 72 vortices (i.e., L = 72), which 
is a half of that for A; = 1 (e.g., when L = 36). In the 
case oi k = 1, a unit cell contains only one vortex. [The 
"unit cell" (UC) is a minimum rcpeatable set of TC(s) 
containing an integer number of vortices. For example, 
for L = 1 the UC is the entire channel with one vortex, 
and for L = 36 the UC is one TC with one vortex.] If 
one vortex can overcome the potential barrier and move 
in the easy direction, all the vortices can do so at the 
same time, i.e., collectively. This results in integer steps 
of the average angular velocity ujq in the to — Iq curve. 
However, fractional steps appear when there are more 
than one vortex in each UC. For instance, a UC contains 
one TC with two vortices when L = 72, and if L = 24, 
the UC is constructed by three TCs with two vortices. In 
general case, the unit cell contains m TCs with k vortices 
inside. If fc > 1, there are more than one vortex in the 
UC, and those vortices are not equivalent, i.e., they are 
not located at the equivalent position in the TC and/or 
they experience different interactions with the boundary. 






FIG. 5: (Color online) Trajectories of vortex motion during one ac period in the case of high density of vortices (L = 80), i.e., 
the multi-vortex regime, for different values of the driving current, Iq = 9 (a), 27 (b), and 36 (c). When the current is in a 
certain range 7o = 5 — 19 [see Fig. [3](c)], the vortices move along 2D trajectory [e.g., Jo = 9, shown in (a)] which is different 
from the circular ID trajectories in the low density case. With increasing current, more vortices start to move along circular 
trajectories and thus 2D trajectories finally turns to a ID circle [e.g., see (c)]. During the transition from the 2D to ID motion, 
the distribution of vortices becomes inhomogeneous [shown in (b)]. 



Therefore, they move with different angular velocities in 
each period. Let us take k = 2, for example (see Fig. 
[5]). As shown in Fig. IH^a), the vortices located in a TC 
that is followed by an empty TC in the easy direction 
(type-A) can overcome the potential barriers in a period 
of alternating current and move to its neighbor TC, while 
the others (type-B) are still localized in its original TC. 
When Iq = 10 [shown in Fig. [SJb)], type-A vortices move 
from C2 to C3 but type-B vortices only oscillate in CI. In 
the view of the whole circular channel, only a half of vor- 
tices (i.e., type-A vortices) move to another TC while the 
other half of vortices (i.e., type-B vortices) do not move. 
After that, the type-B vortices are located at similar po- 
sition as type-A vortices in the previous cycle, i.e., the 
previous type-B vortices now become type-A vortices and 
they will move to next TC in the anticlockwise direction 
during the next period. Therefore, when we calculate the 
average angular velocity for all the vortices, the value is a 
half of the angular velocity of typc-A vortices, i.e., Wo/2. 
If the current increases, every vortex can overcome the 
barriers and move to the neighbor TC [e.g., as shown in 
Fig. [51 type-A vortices move from C2 to C3 and type-B 
vortices move from CI to C2 when Iq = 13]. Then av- 
erage angular velocity becomes wq which is the same as 
that for fc = 1 when all vortices move to the neighbor 
TC. 



B. Commensurability effect of frequency 

Further we analyze the angular velocity evolution while 
varying the ac drive frequency (shown in Fig. |9]). For 
varying frequency v^ the uj — v curves are characterized 
by peaks and/or oscillations. Let us clarify this behavior. 
As discussed in Sec. Ill, if fm > f^ > fin- vortices 



can move in the easy direction and be frozen/blocked in 
the hard direction when the current alternates. Let us 
introduce a time scale, Tg, to characterize the motion of 
a vortex over the entire TC. One TC occupies an angle 
00 ~ 2tt/N . If we assume that a vortex moves over the 
entire TC to the equivalent position in the next TC, then 
To should satisfy the following condition: 



To 



dt = / [f{lQ) +f + r] ■ Odt/r,r = 



(6) 



where 9 is the unit vector in the azimuthal direction. In 
the single vortex case (L = 1), we set f™ = 0, and 
6*0 is a constant. Considering the circular trajectory of 
the moving vortex, the integral of / is also a constant. 
Therefore, Tq only depends on the current Iq and we use 
the notation To(Iq) instead of Tq in order to show the 
dependence on the driving current. 

If /- > /'^(/o) > /+ and the ac period t = 2Tq{Io) = 
to(/o), then the rectified signal will be maximum since 
all the vortices coherently move over the entire TC in 
the easy direction during the first half period and do not 
move backwards in the next half period. Then we can ob- 
tain the principal period io(^o) = 2To(/o) for each value 
of the current and roughly estimate forces /^ and /^ 
from the ui — j^ or cu — t curves [as shown in Fig. [5]. The 
maximum angular velocity Um for different current and 
vorticitics L is shown in the insets of Fig. [S] For example, 
in the case of L = 1 [see Fig. [D^a)], the angular veloc- 
ity is always zero when Iq < 8 and becomes non-zero for 
Iq = 10, i.e., the vortex starts to move in the easy direc- 
tion. Therefore, /'^(/q = 8) < /+ < /''(/q = 10). How- 
ever, for Iq > 22, the maximum angular velocity lj„i de- 
creases as compared to the case for Iq = 20, which means 
/"{Iq = 20) < / - < f^ilQ = 22). For /'^(/q) > /" 
when driving current increases, vortices move backwards 




driving current (I ) 



driving current (I ) 



FIG. 6: (Color online) The CJ — Iq curves for different frequencies: i^ = 0.1 (a), (b), i^ = 1 (c), (d), and i^ = 10 (e), (f). (a), (c), 
and (e) are in the system with a narrower neck, g — 0.1, and (b), (d), and (f) are in the system with a wider neck g = 0.15. The 
scale of Lo in panel (d) is different from the others. For low frequency current (e.g., ;/ = 0.1) the ratchet effect is clearly seen 
[as shown in (a) and (b)]. In this case, the angular velocity first increases continuously to the maximum and then decreases 
to zero. For intermediate frequency v = 1, besides the ratchet effect, the effect of commensurability, i.e., discontinuity in the 
angular velocity uj{Io) (see Sec. IV), has a clear influence on the u — Iq curve [e.g., see (c) and (d) for L = 1 or L = 10]. For 
high frequency current (e.g., v = 10), vortices are confined in their initial TCs and only oscillate in a single potential well. 
Therefore, the effect of the periodic ratchet potential disappears [shown in (e) and (f)]. 



(i.e.; the motion in the hard direction) therefore result- 
ing in the decreasing net angular velocity. For example, 
when /q = 22 [shown in Fig. [H^a)] , the maximum ujm be- 
comes smaller and the jumps in the w — i curves become 
smoother, which means the effect of boundary becomes 
weaker under a strong driving force. Therefore, the an- 
gular velocity should be zero when the driving force goes 
to infinity. It explains why the angular velocity decreases 
and goes to zero when the driving current increases above 
some critical value. 

For the single-vortex case, the frequency dependence 
of the first local maximum ui[^ in the lo — v curves versus 
the corresponding drive frequencies, vq = 1/ioj is plotted 
in Fig. [To] which shows a linear behavior. When the den- 
sity is increased up to one vortex per cell {L/N < 1), the 
a; — t curve is similar to the one for L = 1 for both incom- 
mensurate [e.g., see Fig. IHlJb)] and commensurate cases 



[sec Fig. M,c)] because the interaction between vortices 
is weak. For a higher density, the (in-) commensurability 
effect becomes pronounced. The sharp jumps of the an- 
gular velocity are obtained only in the commensurate 
case [e.g., L = 72 shown in Fig. M,^)] but not in the 
incommensurate case [e.g., i = 80 shown in Fig. IDJf)]. 
Comparing the high density cases with the low density 
cases, we can conclude that sharp jumps in a; — t curve, 
which are both found in single-vortex and multi-vortex 
regimes, are due to the effect of periodically repeated 
boundaries of the ratchet potential and the increasing 
local maximum of the angular velocity when t increases 
in the multi-vortex regime. This is different from the 
single-vortex regime, because of the strong interaction 
between vortices induced by the high vortex density. If 
the principal period (i.e., the time period for the first 
jump) is denoted by to, the others will be harmonics of 
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based on a technique for using weak-pinning channels 
with tailored edges to produce asymmetric vortex con- 
fining potentials. Such an arrangement resulted in sub- 
stantial asymmetric vortex response for linear channels 
on a strip geometry [6|. The Corbino sample consists of 
a Si substrate with a 200-nm thick film of weak-pinning 
a-NbGc and a 50-nm thick film of strong-pinning NbN 
on top. The fabrication followed the scheme of previous 
weak-pinning channel devices p, |30J, with the 1.5-mm 
diameter Corbino disk pattern etched through the entire 
superconducting bilayer. The 500-/im diameter circular 
chain of triangular cells was etched through the NbN 
layer (Fig. [TT|) . thus defining the weak-pinning a-NbGe 
channel region for vortex flow. 



Wirebonds were attached between the center and 
perimeter of the Corbino disk for injecting a bias cur- 
rent with a radial flow. Because of the rather small fiux- 
flow voltages for vortex motion in a single channel, it 
was necessary to use a custom picovoltmeter based on 
a dc SQUID operated in a flux- locked loop [30|- Mea- 
surements of the noise power at different temperatures 
were used to calibrate the value for the series resistance 
at the SQUID input, thus allowing for a measurement of 
the system gain, as described in Ref. [30[. During the 
measurements, the sample and SQUID were immersed in 
a pumped liquid helium bath. Shielding of external mag- 
netic flelds was achieved with a /i-metal shield surround- 
ing the dewar and a superconducting Pb shield around 
the sample and SQUID on the bottom of the cryogenic 
insert. Vortices were introduced into the channel by tem- 
porarily raising the sample above the helium bath, heat- 
ing to -6 K - between T^NbGo ^ 2.88 K and T^^^ « 10 K 
- while applying a small magnetic field with a Hclmholtz 
coil on the insert, then cooling back down below T^^'^'^. 



FIG. 7; (Color online) The ui — Iq curve for a channel with 
g = 0.1 and different k: k ^ 1 (a), 2 (b), 3 (c), and 5 (d). L 
is the vorticity, i.e., the number of vortices. A'^ is the number 
of cells in the circular chain. If L/N = 1/m, e.g., L = 1,36, 
the uj — lo curve shows several steps and the height of each 
integer step is wo- When L/N = k/m (k=2, 3, 5), the height 
of each fractional step is coo/k [e.g., as shown in (b) (c) and 
(d)]. 



to, fcto (fc = 2,3,...). 



V. EXPERIMENTAL DETECTION OF VORTEX 
RATCHET EFFECT IN A CORBINO GEOMETRY 

We have performed preliminary measurements on a de- 
vice containing a single weak-pinning circular channel for 
guiding vortex motion in a Corbino geometry. This is 



Upon reaching the desired measurement temperature, 
the bias current was varied incrementally and the flux- 
flow voltage sensed by the SQUID was recorded for each 
current value. Such current- voltage characteristics (IVC) 
exhibited substantial asymmetries between the critical 
current for vortices to begin to move through the channel 
for the two directions. This asymmetry persisted down to 
the lowest measurement temperature, 1.60 K, well below 
the onset of superconductivity in the channel at T^^Ge 
(Fig. [T^ . The larger critical current corresponded to 
the sense of vortex motion in the hard direction through 
the ratchet cells. Due to experimental wiring limitations 
on these preliminary measurements of a Corbino ratchet 
channel, it was not possible to sweep the bias current with 
an oscillatory drive to study the average flux-flow voltage 
over a cycle. Nonetheless, the large asymmetry between 
the two senses of critical current demonstrates the po- 
tential for a weak-pinning ratchet channel in a Corbino 
geometry to rectify vortex motion. 



FIG. 8: (Color online) (a) The config- 
uration of 24 vortices after 1000 periods 
of oscillations, t = ti. The moving tra- 
jectories of vortices in one period of ac 
current, i.e., from t = ti to t = ti + T, 
is shown for /o = 10 (b), 13 (c), and 15 
(d). The insets of (a)-(d) show the vortex 
distribution in a unit cell. For L = 24, 
a type-A vortex is the vortex in a trian- 
gular cell (TC) [e.g., TC C2 in the in- 
set of (a)] that has a neighbor empty TC 
in the easy direction [e.g., TC C3 in the 
inset of (a)], i.e., a TC without vortex, 
while a type-B vortex is in a TC [e.g., TC 
CI in the inset of (a)] that neighbors a 
TC with type-A vortex inside. When the 
current drives vortices in the easy direc- 
tion (e.g., Jo — 10), two types of vortices 
move in different ways. Type-A vortices 
move from C2 to C3 but type-B vortices 
do not move to C2 due to a larger repul- 
sive interaction in the easy direction (b). 
When Jo = 13 (c), oj = lvq is the same 
as the integer step in a; — Jo curve [see 
Fig. Ob)]. The type-A (type-B) vortex 
moves over one TC and remains type-A 
(type-B) in the end of the period. If the 
current increases further, e.g., Jo = 15, 
type-A and type-B vortices move with 
different net angular velocities and af- 
ter every period they switch their type 
[shown in (d)]. 




VI. CONCLUSIONS 

A vortex moving in an asymmetric eircular channel 
in a Corbino setup experiences the confinement of the 
boundary, the repulsive interaction due to other vortices 
and the gradient (i.e., radially decreasing) driving force 
when an external current is applied. The combination of 
these factors determines the vortex motion. Different dy- 
namical behavior is observed in low and high density sys- 
tems, which are referred to as "singlc-vortcx" and "multi- 
vortex" regimes, respectively. For low density, i.e., in the 
single- vortex regime, there is no more than one vortex 
per triangular cell so that the vortex-vortex interaction 
is negligible. Therefore, the ratchet potential due to the 
boundary dominates and all the vortices follow circular 
ID trajectories. Since all the vortices still move in ID 
when there is more than one vortex per cell (i.e., a higher 
density) , a vortex can escape even easier either in the easy 
or in the hard direction due to the interaction between 
vortices in the same cell. This results in decrease of the 
rectified net current. However, when the number of vor- 
tices increases further (i.e., in the "multi- vortex" regime 
of rectification when the rectified net current increases 
with the vortex density), the trajectories of vortex mo- 
tion for low driving currents are 2D while these trajecto- 
ries squeeze and turn to ID circle with increasing driving 



current. Because of the circular geometry of the chan- 
nel, the density of vortices first becomes inhomogeneous 
[i.e, see Fig. ^h)] during the transition from 2D motion 
to ID motion in the multi-vortex regime and then be- 
comes homogeneous again when all the vortices move in 
a circular trajectory. 

Considering the asymmetry in the radial direction, vor- 
tices near the inner/outer corner of the triangular cells 
(TCs) (i.e., closer/further to the center of the disk) are 
driven by different Lorentz forces and for some specific 
value of driving current the vortex in the outer corner 
moves to the inner corner while the one in the inner cor- 
ner moves to the next TC [e.g., see Fig. EJa)]. This kind 
of motion prevents vortices from arriving simultaneously 
at the narrow part that would lead to jamming which 
occurs in case of a linear channel. When the density in- 
creases, the maximum net angular velocity ujm remains 
the same in the single- vortex regime and then decreases 
until reaching the multi- vortex regime, and then w„j in- 
creases. 

The frequency of driving current also strongly infiu- 
ences the vortex dynamical behavior. The ac frequency 
determines the possible distance a vortex moves during 
an ac period. For high frequency, a nearly zero net an- 
gular velocity is obtained for different values of driving 
current Iq . Each vortex is unable to overcome the energy 
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FIG. 9: (Color online) The u; — t curves in a channel with g = 0.1 for different density of vortices: L/N ^ 1 (a), L/N < 1 (b), 
L/N = 1 (c) (commensurate), 2 > L/N > 1 (d), L/N = 2 (e) (commensurate), and L/N > 2 (f). When the density L/N < 1, 
i.e., in the single-vortex regime, the angular velocity uj{t) first reaches the maximum and then oscillates [see (a), (b) and (c)]. 
With increasing current, the maximum velocity increases until the driving force reaches the value larger than /„. For high 
density, the angular velocity ij{t) does not reach the maximum because the vortex motion in the easy direction is compensated 
by that in the hard direction. uj{t) increases when the frequency decreases. In commensurate cases [(c) and (e)] a;(t) oscillates 
with sharp jumps which are not observed in incommensurate cases for high vortex density (f). 



barrier and is confined by a single potential well instead 
of the periodic ratchet potential. When the frequency is 
low, the ratchet effect is clearly observed in the w — /q 
curve but the commensur ability effect is not present. For 
an intermediate frequency of driving current under which 
the distance of moving vortex in a period is comparable 
to the size of the TC, both the ratchet effect and the 
commensurability effect have been observed. 

Besides the vortex density and the frequency of cur- 
rent, the commensurability between the numbers of vor- 



tices and TCs also plays an important role in the dy- 
namical behavior of vortices, which leads to jumps in the 
angular velocity w with increasing driving current Iq (i.e., 
steps in the uj — Iq curve). Therefore, under some specific 
conditions, the average angular velocity of vortices is not 
a continuous function of the driving current. The com- 
mensurability also influences the minimum difference of 
the angular velocity for different steps in the uj~Iq curve, 
and results in integer steps (i.e., the large steps which are 
found, e.g., for vorticity L = 1) and fractional steps (i.e.. 
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the smaller steps whose magnitude are fractions of the 
magnitude of the integer steps, e.g., for vorticity L = 24) 
in a certain range of the current. We also obtained several 
peaks (sharp jumps) in the uj — t curve, which correspond 
to the principal period Iq (during which the vortex can 
move over one TC in the easy direction), and harmon- 
ics periods /cj/q {k — 2, 3, ...) (during which vortices move 
over k TCs in the easy direction). The net flow of vortices 
is enhanced when the ac period is one of the harmonics 
periods, i.e., the average angular velocity reaches a local 



maximum. 



FIG. 10: (Color online) The frequency dependence of the first 
local maximum uj'^-^ of the u)~t curves for vortices in the single- 
vortex regime when the driving force is smaller than /^ . The 
magnitude of the first peak shows a linear dependence on the 
principal frequency vq. 



*' 4.5 \mr^ 




FIG. 11: (Color online) (a) Detail of single ratchet cell of 
channel and schematic of Corbino ratchet chip layout; extra 
pads and leads beyond disk were not used for the measure- 
ments presented here, (b) Atomic force microscope image 
of portion of Corbino ratchet channel, (c) Scanning electron 
micrograph of portion of Corbino ratchet channel. 
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FIG. 12: (Color online) Current-voltage characteristic of 
Corbino ratchet channel cooled to 1.60 K in an external mag- 
netic field of 0.26 Oe as described in text. 
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